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, Abstract. Standard GRB fireballs must carry free neutrons. This 

' crucially changes the mechanism of fireball deceleration by an external 

■ medium. As the ion fireball decelerates, the coasting neutrons form a 

<D ' leading front. They gradually decay, leaving behind a relativistic trail of 

decay products mixed with the ambient medium. The ion fireball sweeps 
up the trail and drives a shock wave in it. Thus, observed afterglow emis- 
sion is produced in the neutron trail. The impact of neutrons turns off 
at ~ 10 17 cm from the explosion center, and here a spectral transition is 
expected in GRB afterglows. Absence of neutron signatures would point 
l/-) ' to absence of baryons and a dominant Poynting flux in the fireballs. 
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O . 1. Introduction 
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Importance of neutrons in GRBs was realized recently (Derishev, Kocharovsky, 
& Kocharovsky 1999a,b; Bahcall & Meszaros 2000; Meszaros & Rees 2000; 
Fuller, Pruet, & Abazajian 2000; Pruet & Dalai 2002; Beloborodov 2003a,b). 
Their presence is inevitable in a standard fireball (sections 2 and 3) and they 
qualitatively change the explosion picture at radii of 10 16 — 10 17 cm (section 4). 
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^ . 2. Neutronization of the Central Engine 

h : 

GRB central engines are very dense and hot (e.g. Meszaros 2002). The central 
matter (or matter surrounding the central black hole) is filled with hot e^ pairs 
in thermodynamic equilibrium with blackbody radiation at a temperature kT = 
1 — 10 MeV. Its baryonic component is made of neutrons and protons (nuclei 
break up into free nucleons in the unshadowed region of the T — p plane shown 
in left panel of Fig. 1) and frequent e^ captures on nucleons take place, 

e~ + p — >n + v, e + + n — > p + v. (1) 

These reactions quickly convert protons to neutrons and back, and establish an 
equilibrium proton fraction Y e = n p /(n n +n p ) which is shown in Figure 1. Any 
plausible GRB engines belong to the T — p region where Y e < 0.5. So, the 
neutrons are not only present in the central engines — in fact, they dominate. 

When the neutronized matter is ejected into a high-entropy fireball, it re- 
mains neutron-dominated because (1) expansion occurs too fast and Y e freezes 
out at its initial value below 0.5, or (2) the fireball partially absorbs the neutrino 
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Figure 1. Contours of equilibrium Y e = n p /(n n + n p ) on the T-p 
plane for z^-transparent (left) and ^-opaque (right) matter. Thick curve 
shows the boundary of free-nucleon region. (From Beloborodov 2003a) 

flux from the engine, which tends to keep Y e < 0.5 (see Beloborodov 2003a for 
details). 

3. Ejected Neutrons Survive the Nucleosynthesis 

As the fireball expands and cools, nucleons tend to recombine into a-particles 
(like they do in the big bang). This process competes, however, with rapid ex- 
pansion and can freeze out. For this reason, nucleosynthesis is suppressed in 
fireballs with a high photon-to-baryon ratio (f) = n^/rib (or, equivalently, high 
entropy per baryon s/k = 3.60). A minimum 4> in GRBs is ~ 10 5 , which is just 
marginal for nucleosynthesis (detailed calculations of nucleosynthesis in GRBs 
have been done recently by Lemoine 2002; Pruet, Guiles, & Fuller 2002; Be- 
loborodov 2003a) . Even in the extreme case of complete nucleon recombination 
there are still leftover neutrons because of the neutron excess (Y e < 0.5). Their 
minimum mass fraction is X n = 1 — 2Y e . 

In addition, synthesized helium is likely destroyed during the subsequent 
evolution of the explosion. This happens if (1) there appears a substantial 
relative bulk velocity between the neutron and ion components (as a result of 
neutron decoupling during the acceleration stage of the fireball) or (2) internal 
shocks occur and heat the ions to a high temperature (Beloborodov 2003a). 

Thus, a substantial neutron component appears inevitably in the standard 
fireball scenario. Neutrons develop a Lorentz factor T n = 10 2 — 10 3 at the very 
beginning of the explosion when the fireball is accelerated by radiation pressure: 
they are collisionally coupled to the ions in the early dense fireball, and decouple 
close to the end of the acceleration stage. Then the neutrons coast and gradually 
decay with a mean lifetime Tp ~ 900 s and a mean decay radius Rp = crpT n , 
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4. Neutron-Fed Blast Wave 

Let us remind what happens in a relativistic explosion without neutrons. The 
ejected fireball with mass M e j and Lorentz factor T c j sweeps up an ambient 
medium with density uq = 1 — 100 cm -3 and gradually dissipates its kinetic 
energy. The dissipation rate peaks at a characteristic "deceleration" radius 
Rd 

ec ~ 10 cm where half of the initial energy is dissipated. Rdec corresponds 
to swept-up mass m& ec = M e j/r e j. Further dynamics is described by the self- 
similar blast wave model of Blandford & McKee (1976). How does this picture 
change in the presence of neutrons? 

At radii under consideration, R > 10 15 cm, the ejected fireball is a shell of 
thickness A <C R. In contrast to neutrons, the ion component of the fireball is 
aware of the external medium and its Lorentz factor T decreases. As T decreases 
below T n , the ions fall behind and separate from the neutrons. Thus the fireball 
splits into two relativistic shells which we name N (neutrons) and I (ions). The 
mass of the leading N-shell is decreasing because of the /3-decay, 

M n (R)=M^exp(--^j. (3) 

The N-shell energy, E n = T n M n c 2 , is huge compared to the ambient rest- 
mass mc 2 even at R > Rp. For example, at R = R^ cc we find E n /m,i ec c 2 = 
A n r n r e j exp(— Rfcc/Rp) where X n = M®/M e - S is the initial neutron fraction of 
the fireball. 

The neutron decay products p and e~ share immediately their huge momen- 
tum with ambient particles due to two-stream instability and the N-shell leaves 
behind a mixed trail with a relativistic bulk velocity (3 < (3 n (Beloborodov 2003b) 

m = 1 + atc)- 1 ' 7(jR) = (i-p 2 ) 1 / 2 = W+^c+W (4) 

where 

C{R) = ^m~ = R^ [dRj ' (5) 

and m{R) is ambient mass enclosed by radius R. There exists a characteristic 
radius .Rtraii where the trail becomes nonrelativistic (f3 = 0.5). It is defined by 
condition ( = T~ l , which requires about 10 e-folds of the decay (for a typical 
mp ~ mdcc ~ 1(T 5 M0). Thus, 



#trail ~ 10^ = 0.8 x 10 1 ' cm. (6) 

^trail depends very weakly (logarithmically) on the ambient density and the 
initial neutron fraction of the fireball, X n . 

The decaying N-shell not only accelerates the medium as it passes through 
it. It also compresses the medium, loads with new particles, and heat to a high 
temperature. The rest-frame density and relativistic enthalpy of the trail are 

n = „ (i + < )(^ + 2 r„ C + i) 1/2 , + + (7) 



4 



Beloborodov 



where uq is the medium density ahead of the N-shell. For T" 1 < ( < T n one 
finds jU S> 1, i.e. the thermal energy of the trail far exceeds its rest-mass energy. 

The ion fireball follows the neutron front and collects the trail. As a result, 
(1) the ion Lorentz factor T decreases and (2) a shock wave propagates in the trail 
material. The shock has a Lorentz factor T s h ^ T and it cannot catch up with the 
neutron front (unless rio[R] falls off steeper than R~ 3 ). Dynamics and dissipation 
in the shock are discussed in Beloborodov (2003b). We emphasize here important 
differences from a customary external shock: the neutron-trail shock propagates 
in a relativistically moving, dense, hot, and possibly magnetized medium behind 
the leading neutron front. 

The neutron impact ceases at -Rtraii ~ 10 17 cm, which can leave an imprint 
on the observed afterglow. For example, the shock dissipation can have a second 
bump (Beloborodov 2003b), and a spectral transition is also possible. The arrival 
time of radiation emitted at i?traii is approximately i?traii/2r 2 c (counted from the 
arrival of first 7-rays). It may be as long as 30 days or as short as a few seconds, 
depending on the fireball Lorentz factor r(i? trai i). Recent early observation of a 
GRB afterglow (GRB 021004) discovered an interesting re-brightening at 10 3 s. 
Also, we do not exclude a possible relevance of neutrons to the 20 day bumps 
observed in a few GRBs, as the time coincides with -Rtraii /c. 

Neutron signatures should be absent if the fireball is dominated by a Poynt- 
ing flux and has extremely low baryon loading. Then the neutron component 
decouples early, with a modest Lorentz factor T n , and decays at small radii. The 
upper bound on T n due to decoupling is T n « 300(M n /10 26 ) 1 / 3 where M n [g/s] 
is the mass outflow rate per unit solid angle of the fireball (Beloborodov 2003a). 

An additional piece of GRB physics is the interaction of the prompt 7-ray 
radiation with an ambient medium (see Beloborodov 2002 and refs. therein). It 
leads to a gap opening and crucially affects the early afterglow emitted at radii 
R < 10 16 (£' 7 /10 53 ) 1 / 2 cm, where E 1 [erg] is the isotropic energy of the GRB. 
Here, we focused on larger radii where the neutron effects dominate. 
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